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Abstract. The Hough transformation has been used suc- 
cessfully for more than four decades. Originally used for 
tracking particle traces in bubble chamber images, this 
work shows a novel approach turning the initial idea into 
a powerful tool to incoherently detect millisecond pul- 
sars in binary orbits. 

This poster presents the method used, a discussion on 
how to treat the time domain data from radio receivers 
and create the input "image" for the Hough transforma- 
tion, details about the advantages and disadvantages of 
this approach, and finally some results from pulsars in 47 
Tucanae. 
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Introduction 

The problem of discovering millisecond pulsars in binary 
orbits is apparent. Not only the parameter range for iso- 
lated pulsars has to be addressed, but — at least — five 
Keplerian parameters have to be included additionally. A 
simple Fourier analysis is not sufficient for this class of pul- 
sars, since the signal is smeared out over a large frequency 
band due to Doppler's effect (ref. Figure. EJ). 
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Fig. 1. Power spectrum of the binary pulsar 47 Tuc J. 



Over the past decades quite a few search algorithms 
have been devised to counter this problem, e.g. the very 
successful "acceleration search" [T]or the sideband searchO 
However, these approaches work purely on coherent data 
sets, thus they are ultimately limited by the maximum 
possible length of an observation 

We would like to present an alternative approach to 
search for weak pulsars in a binary orbit, which utilizes 
data taken in several observation runs. In the following 
section we will briefly explain how the Hough transforma- 
tion works, and how it can be applied to search for binary, 
millisecond pulsars before summarizing some of the results 
obtained for 47 Tucanae. 



The Hough transformation 

Paul Hough developed and patented the algorithm around 
1960. The original algorithm was designed to search for 
straight particle traces in bubble chamber images by sim- 
ply turning the usual line equation y = ax + b into a 
parameter space equation b = y — ax. 

From these two equations it is obvious that a line in 
the original space ("xy-space") corresponds to a single 
point in the parameter space ("ab-space"). On the other 
hand a single point in xy-space is represented by a line 
in ab-space. Although sounding trivial, it actually solves 
a complex problem: Identifying parameters of lines found 
in (pixel) images. 

Figure [2] shows what happens when transforming many 
points from a single line. For each of these points a line is 
drawn in ab-space and since each of these points belong 
to the same xy-line, their lines are required to intersect 
in a single point in ab-space. Naturally, random points in 
xy-space convert to random lines in ab-space only adding 
up to a random background. 

Tracking millions of lines along with their intersections 
is costly in terms of computing power and memory, thus 
we chose to discretize the ab-space into rectangular pixels 
and simply count how many lines are drawn through each 
pixel. This approach allows us to discover the blue "line" 
in Figure[3l where only 95 signal points are hidden beneath 
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Fig. 2. Points from a line in xy-space transform into lines 
intersecting in a single ab-point. 

200,000 noise points while still be detectable at an SNR 
of about 20! 

For binary pulsars, the frequency signature can be de- 
scribed by 

/(t)=/p(l + A COs(27T/ -* )), 

where f p is the intrinsic frequency of the pulsar, f Q the 
orbital frequency, A the (dimensionless) frequency ampli- 
tude due to Doppler modulation and <£> the orbital phase 
at t = 0. This equation can be recast into a line equation 

C = Ax + By, where 

A = sin(27r/ -$ ), 

B = cos(2tt/ - $ ), 

C = (f(t) - f p )/f p , 

x = / p A sin<l> , 

V = fp^o cos$ . 

Here, we split the variables into two types: 

A , 3> : These variables are considered in the Hough 
transform, i.e. the resulting Hough planes (formerly 
ab-space) are polar diagrams in these parameters. 
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Fig. 3. Recovering a hidden line (top, blue) with a SNR of 
about 20 (tiny white point (bottom)). 



f P , fo' Called "external" variables, because they appear 
only as input parameters to the Hough transform and 
are varied outside the scope of it. This separation al- 
lows us to distribute the work onto a cluster of com- 
puters. 

So far, we have not said anything about our "input 
image" . From the time domain data0 we use small pieces 
of the data (typically a few minutes), create a power spec- 
trum for each piece and select peaks from the spectra. 
Then, we draw the registered peaks into a time-frequency 
plane which serves as our input for the Hough transfor- 
mation. 



47 Tucanae results 

In this section we summarize some results obtained from 
47 Tucanae data. Figure [H shows the outcome obtained 
for the full parameter space (f p , P Q = f~ l ). 

The maximum number counts (of how many lines are 
drawn through any pixel) found in each of the about 10 9 
computed Hough planes are color coded in this image. 
One can clearly see the signatures of several noise sources 
(lines at and below 100 Hz as well as pulsars (173 Hz, 
186 Hz) and their harmonics at twice, thrice, ... of the 
fundamental frequency. 

A closer look at the Hough planes reveals that 47 Tu- 
canae C at 173.709 Hz is an isolated pulsar (Figure [5}. 
Please recall that such a source does not feature any no- 
ticeable Doppler modulation in the solar system barycen- 
ter and the parameter A is therefore equal to zero. This 
in turn is shown by the maximum at the center/origin. 

For the binary pulsar 47 Tucanae J the situation is 
different. It orbits its companion in less than 3 hours and 



1 Kindly made available for us by the pulsar group at Jodrell 
Bank, Manchester, UK 
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Fig. 4. 47 Tucanae global result map 



173.709Hz, 101052 s, max at 0.26 mHz and 7.72 deg 




Fig. 5. 47 Tucanae C: The maximum in this Hough plane 
is at the origin, suggesting no Doppler modulation of the 
signal (isolated pulsar). 



the maximal measurable Doppler shift for this system is 
about 12 mHz. Figure [6] shows the Hough plane featuring 
the maximum at a distance to the origin corresponding 
exactly to this Doppler shift. 

More already known pulsars have been found in the 
data set, but are not displayed here, please refer to [3] for 
more information. 



Summary 

Although this new method has not revealed any new pul- 
sars in the well observer globular cluster 47 Tucanae (yet), 
many already known pulsars have been "rediscovered" in 
a blind search. Since the obtained parameters for these 
sources agree remarkably well with the already published 
values, the method has proven to be principally able to 
discover new pulsars. 
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Fig. 6. 47 Tucanae J: The maximum in this Hough plane 
is not at the origin anymore. The offset indicates a Doppler 
modulation of about 12 mHz agreeing well with literature 
values. 

The most important advantage over "classical" search 
algorithms is the semi-coherent approach, allowing to 
combine data taken on several occasions into a single, 
more sensitive search. However, it should be noted that 
the Hough transform requires a large amount of CPU 
power, a typical search over a frequency bandwidth of a 
few hundred Hz can easily take several 10,000 CPU-hours 
- the current focus is to lower this number. 

Finally, the author would like to thank the pulsar 
group at the Jodrell Bank observatory, especially Michael 
Kramer and Dune Lorimer, and the Albert-Einstein- 
Institute for the opportunity to work on this project over 
the past years. At AEI especially the members of the data 
analysis groufH proved to be invaluable sources of knowl- 
edge. Special thanks go to Maria Alessandra Papa, Curt 
Cutler, Bernard Schutz, Alicia Sintes, Steffen Grunewald, 
Bernd Machenschalk, Badri Krishnan, Reinhard Prix and 
many more! 
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